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ABSTRACT 

We investigate the large behavior of the complex solutions for the two magnon 
system in the S=l/2 Heisenberg XXZ model. The Bethe ansatz equations are nu- 
merically solved for the string solutions with new iteration method. A clear evidence 
of the violation of the string configurations is found at = 22, 62, 121, 200, 299, 417, 
but the broken states are still Bethe states. The number of the Bethe states is con- 
sistent with the exact diagonalization, except one singular state. 
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1 Introduction 



In 1931, Bethe presented a unique and genius way of solving one- dimensional Heisenberg 
model with spin S = 1/2 systems |]l], The structure of the Bethe ansatz solutions has 
been studied quite extensively. Since the Bethe ansatz method has a history of many 
years, the fundamental properties of the Bethe ansatz solution for the Heisenberg XXZ 
model are well understood by now p, ^ |], |^. However, the completeness of the Bethe 
ansatz solutions is still discussed actively p|, ^, |l^ . Bethe ansatz solutions contain complex 
solutions whose rapidities are complex values. For these solutions, Bethe proposed the 
string hypothesis. He assumed the shape of the complex solutions which are called " string 
solutions" and showed that the real solutions together with such complex solutions give 
the correct number of the states [|I| . Therefore, it is widely believed that the completeness 
of the Bethe's state holds true. 

In two magnon system, the solutions are parametrized by two integer numbers mi 
and m2- The string solutions can be classified into two different types, depending on the 
parity of mi + m2 in the two magnon state. In fact, the behaviors of the wave function for 
the two string solutions are quite different from each other. The string solutions with odd 
parity of mi + m2 were first discussed by Bethe, and since then they have been discussed 
quite extensively. In particular, Essler, Korepin and Schoutens (EKS) carried out the 
careful study of the string solutions. Here, we call this type of the string "EKS-string" 
as we define later in detail. Recently, it is found that this EKS-string solutions break 



down to two real solutions at = 22,62, 121... ||12|, |T3|, |14|, which are partly predicted 
by Essler, Korepin and Schoutens. On the other hand, there are other type of complex 



solutions which are treated by Vladimirov [Tl|. This string solutions have even parity of 
mi + m2, but their behaviors are quite different from the EKS-strings. Here, we call them 
"V-string". According to the string hypothesis, the imaginary part of the string is 1/2. 
However, Vladimirov shows that the imaginary part of the V-strings behaves like 

In this paper, we present the large behavior of the string solutions of the Bethe 
ansatz equations for the two magnon system. Here, we have solved the Bethe ansatz 
equations numerically up to quite large number of the site A^. The numerical evaluation 
has some difficulties for the complex solutions in the Bethe ansatz equations. A simple 
iteration method cannot give any convergent results, and thus we have to develope a new 
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way to solve them. This is what we have achieved here as we discuss it later in detail. The 
solutions are compared with those of the exact diagonalization method up to = 180. 
Here, we indeed confirm the violation of the EKS-strings up to = 417. Further, they 
are not out of the Bethe Ansatz solutions, and therefore, the number of the state is 
unchanged. In addition, we confirm that the behavior of the V-string rapidity predicted 



by the 1/A^ expansion method [TT[] is indeed consistent with our numerical solutions. 

The paper is organized as follows. In the next section, we briefly explain the solutions 
of the Bethe ansatz equations. In section 3, we present a new iteration method of solving 
the Bethe equations numerically. In section 4, the string solutions for large A^ cases are 
discussed and compared with the predictions of the 1/A^ expansion method. Section 5 
summarizes what we have clarified. 



2 The solutions of the Bethe Ansatz equations 

The Heisenberg XXZ model with spin 1/2 is a model which is most frequently studied 
among all the models of the spin systems. It is solved by the Bethe ansatz technique. 
Here, we briefly describe the Heisenberg model and the Bethe ansatz solutions. The XXZ 
model is described by the following Hamiltonian 

^ . J - - 

H = J (^SfS^_^_l + SfSf^i + AS^S^^i) = — [S^Si^i + + 2AS'j S^^i) , (1) 

1=1 i 

where S'f is a spin operator at the site i and S^{S^) flips down (up) the spin, and A is 
the anisotropic parameter. The periodicity Sn+z = Si is assumed. For XXX model, we 
take A = 1. This Hamiltonian can be diagonalized by the following Bethe state for 
two magnon system, 

\^2)= M^uX2)S-^S-^\0), (2) 

Xl<X2 
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where |0) is the ferromagnetic state with all spins up. The coefficient A{xi, X2) is assumed 
to be of the following shape, 

and must satisfy 

A{xj, Xj) + A{xj + 1, Xj + 1) - 2A{xj, Xj + 1) = 0. (4) 

Therefore, the phase shift (p should satisfy the following equation, 

. . ki- h 
if _ Asm^— 

COS — A cos — - — 

2 2 

where (p is taken to be — tt < (fi < tt. Further, imposing the periodic boundary conditions, 
we have 

kiN + ip =27rmi, (6) 
k2N — (p — 27rm2, (7) 

where nii and 7712 are integers running between and — 1. Without loss of generality, 
we can take mi < 1712- In this case, the energy eigenvalue of the Hamiltonian can be 
written as 



2 

E/J=J2coskj + A(^-2Y 
.1=1 ^ ^ 



(8) 



3 

We can define the rapidity as 

1 k 



A, = -^cot^. (9) 



In this case, the energy is given by 



7 = 1 ■> ' ^ ' 
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In order to obtain the energy eigenvalues of the Hamiltonian, one has to solve eqs.(^), (^) 
and (^. Here, we show the rapidity and the energy in the case of = 8 and A = 1 for 
simplicity. We must solve the equation 

COS -^{'mi + 1712) COS I — 77(,"^i "^2) = COS I — - — K — —{mi — m2) , (11) 



where K = {ki — k2)/2. This equation can be solved analytically once we specify the values 
of mi and m2. It is difficult to prove the completeness of the Bethe Ansatz solutions. 
Therefore, we must carefully treat the Bethe Ansatz equation whose solutions are indeed 
the answer of the model. 

In table |I], we show the calculated results of the energies for A^ = 8 with the two down 
spins together with the energies by the exact diagonalization for A = 1. The detailed 
method to solve the Bethe equations numerically will be given in section 3. 

As can be seen from table 0, there is one state, denoted by *, which cannot be 
reproduced by the Bethe ansatz. The configurations of the states are different for different 
values of A. Here, we only show the results for A = 1 case. For other values of A, we 
only make a comment here. We have the real solution for A = 1 in the Category I, which 
is defined in 0. The string solutions appear for the Category II and HI. For A = 2, there 
is an irregular solution at (mi,m2) = (0,7). Even though the element belongs to the 
Category I, the solution is a string. We have a similar configuration for the case A > 1. 
On the other hand, there is no complex solution for A = 1/2 at A^ = 8. In what follows, 
we will only discuss for A = 1 case. The general A cases will be treated elsewhere. 

Now, there is a state which cannot be reproduced by the Bethe ansatz solution for 
A = 1. We can easily verify that the energy of this state is given by 

E/J = A (^^ - l) . (12) 

The corresponding state is 

N 

\Non - Bethe) = J](-)^5rSr^^ |0) . (13) 

i=l 

Note that the kinetic term of this state is null, 

J2 {S^Sr^i + SrSti) \Non - Bethe) = 0. (14) 
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Thus, we cannot describe this state in the Bethe Ansatz equation because eq.(^ is not 
vahd for this case any more. We note here that this state is not reproduced by the 
algebraic Bethe Ansantz at the isotropic point either WE . 



Bethe suggested that there is a string solution whose imaginary part is infinity 
However, the wave function is divergent. For appropriate normalization, we get the wave 
function given as eq. ([T3|) . Therefore, the wave function eq.([T3|) corresponds to the singular 
Bethe state. 

We comment on the degeneracy of the energy states. Suppose that the ki,k2 are 
solutions for the configuration (mi, 7712) (1 < mi < 1112 < N — 1) with the phase ip. The 
Bethe equations are 

kiN + ip = 2-nmi, k2N - (p = 2Tcm2. (15) 

Now we consider the configuration (A^ — m2, A^ — mi) with k'l, k2,ip'. We can easily show 
that two states have the same energy if the following equations are satisfied 

k[ + = Att - {ki + k2), k[ - k'^ = ki - k2, 'p' = 'P- (16) 

Note that 

= ki {i = 1, 2), for mi + m2 = A^. (17) 

In this case, there is no pair and the Bethe state is a single state. For (0,m2), (1 < m2 < 
A^ — 1) case, the configuration (0, A^ — m2) is a solution with the same energy as the 
configuration (0, m2) when 

k[ + k'2 = 27c-{ki + k2), k[- k'2 = -{ki- k2) -27T, <p' + <p = 0. (18) 



3 New iteration method 

The Bethe ansatz equations (5), (6) and (7) can be solved by the iteration method even 
for the large A^ cases. The convergence of the iteration method is good for the normal 
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Table 1: The energies E^xact and Esethe which are solved by the exact diagonalization 
and by the Bethe ansatz method are presented for the N — S and A = 1 case. * means 
that there is no corresponding state in the Bethe Ansatz method. 



Eexact/J -E'Bethe/'^ 



(mi,m2) 



fe2 



-1.801938 -2cos- 
7 



(3,5) 



67r/7 



87r/7 



67r/7 



-1.267035 



-1.144123 



0.2928932 
0.4370160 

1.0 

1.0 

1.144123 

1.246980 
1.525687 

1.707107 



-1.267035 



i/IO-f t/2 



-0.4450419 -2cos— TT 
7 



10- 



-0.4370160 



-0.2586520 -0.2586520 



2 - V2 



%/lO - V2 



* 



5 

-2 cos -TT 

7 



1.525687 



2-1- V2 



2,4- 
4,6 

2,5 
3,6 

2,6 

1,4 
4,7 

1,5 

3,r 

0,4 
1,3 

5,7' 

0,3 
0,5 

1,6 

2,7' 

0,2 
0,6 



1,2 
6,7; 

1,7- 

1,1 
7,7 

0,1 
0,7' 



0.60357r 
1.10357r 

0.58757r 
0.83757r 

47r/7 

0.31847r 
1.06847r 

0.30857r 
0.80857r 



7r/3 
47r/3 




0.29907r 
0.54907r 






0.89657r 
1.39657r 

1.16257r 
1.41257r 

107r/7 

0.93167r 
1.68167r 

1.19157r 
1.69157r 

7r 

27r/3 
57r/3 

37r/4 
57r/4 

1.45107r 
1.70107r 

7r/2 
37r/2 



0.82797r 
0.82797r 

0.70037r 
0.70037r 

47r/7 

0.54757r 
0.54757r 

0.46847r 
0.46847r 



27r/3 
27r/3 




0.39207r 
0.39207r 






37r/8 -I- 0.9578i 37r/8 - 0.9578i tt -|- 7.6626i 
137r/8 -I- 0.9578i 137r/8 - 0.9578i n + 7.6626i 



27r/7 



127r/7 



7r/4-|-0.3945i 7r/4 - 0.3945j 
77r/4 -I- 0.3945i 77r/4 - 0.3945i 



7r/4 
77r/4 



27r/7 

3.1559i 
3.1559i 






2.0 



(0,0) 
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cases, but the complex solutions cannot be obtained easily by the simple-minded iteration 
method. Here, we present somewhat a different way of numerically solving the Bethe 
equations. But before going to the new method, we briefly explain the usual iteration 
method. The (p in eq.(5) can be written as 

A sin ^2—^ 

COS A cos 

2 2 

The usual iteration method is done as follows 

= (20) 

where wc start from some initial value of Lp^^\ k\ and are given in eqs.(6) and (7). As 
mentioned above, this procedure works well for the normal cases. But for the complex 
solutions, it does not work. 

In this paper, we show a new way of solving the complex solutions. First, we define a 
new variable v by 

ip = iv for {m,m) (V-string), (21) 

ip = t: + iv for (m, m + 1) (EKS-string), (22) 

where v is taken to be real. For V-string case (m, m), we obtain 
coth"- -^^'""^ 



2 2rmT v 

cos A cosh — 

N N 



2m7r . , V 

sinh - = — ^ — cosh - = X{v)^ 

2 —A sinh — 2 

v^2\n(x + VA:2 + 1^ . 
Therefore, the iteration equation becomes 

v^^ = 2 In f X + VA:2 + l) . (23) 

\ / t;=j;(»-l) 



8 



On the other hand, for the case of EKS-string (m, m + 1), we have 



^ -Asmli — 

tanh - — 



2 (2m+ l)7r ^ v ' 

cos — A cosh — 

N N 



(2m + 1)7T . , V 

V Acosh— ^ 

cosh - = — — sinh - = Y(v), 

^ —A smh — ^ 

N 

v^2\n(Y + VY^ - . 

Therefore, the iteration equation becomes 

v^^ ^ 2\n (y + VY^ - l) . (24) 

We can carry out the iteration procedure for miUion times, and obtain good convergent 
results. In this way, this iteration method gives complete solutions, and thus we can find 
all of the string solutions of V-string as well as EKS-string. With this method, we can 
find the string solutions up to ~ 5000. 



4 String solutions for large A^^ 

Next, we discuss the properties of the string solutions. For the string hypothesis, the 
complex solutions are given for M-magnon system by 

Xj^x + i (^^^ - + 0(e-^^), (j = 1, 2, • • • M), (25) 

where x is the same real part of the complex roots \j, and 5 is a positive parameter. 
According to the string hypothesis, the imaginary part y of the string at large N is given 
as, 

y - 1 (26) 
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for the two magnon case (M = 2). On the other hand, Vladimirov predicts |Tl|] the 
non-string type complex solutions by making use of the expansion method. Defining 
the real and imaginary parts of the rapidity hj X = x ±iy, they can be expressed as 



_ iV / 1 TlH^ 



+ 



7l£ 



1 



~2A 



where i = 2,4,6, ■■■< y N. For large A^, they behave as 

X ^ N, y a/zV. 



(27) 



(28) 



In table H, we show the complex solutions of the XXX model for several cases of N. The 
Vladimirov's solutions are given in the following configurations, 



(mi, 7712) = (r, r), r 



1,2,--- 



(29) 



where [x] denotes the maximum integer n for n < x (Gauss's symbol). In figures |1] and 
1^, we show the large behavior of the imaginary part and the real part of the complex 
solutions, respectively. We can see that the large N behavior of the configuration for 
(mi, 7712) = (1, 1), (2,2) is indeed consistent with eq. (pHj). 

On the other hand, the large N behavior of the configurations for m2 = mi + 1 case is 
different. We show the behavior of the imaginary part and the real part of the complex 
solutions in figure ^ and figure There are branch points where the string solutions 
break down. At the branch point, the imaginary part is close to zero, and the real part 
is divided into two parts. However, the broken string becomes two real solutions in the 
Bethe Ansatz equation and thus these broken states are still degenerate, which is shown 
in table |^. In table ^, we list the configurations for the complex solutions. 

In figure I, the minimum number of mi and the sites number A^ are given. The number 
of the string does not increase linealy for A^. Note that, in both cases, there are string 
solutions whose imaginary part is 1/2. However, the value of the imaginary part starts 
to deviate from 1/2 when the value of the site number A^ increases, and finally becomes 
zero for sufficiently large values of A^. 

Finally, the string configurations for the XXX model (A = 1) are classified by the 



large A^ behavior. In |T^, Ilakovac et al. distinguish the string solution in terms of the 
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Table 2: The energies E^xact and the string solutions and non-string type solutions which 
are solved numerically for the XXX model (A = 1). X^iad is the Vladimirov's prediction. 

N Eexact/J (mi,m2) A Klad 



12 2.76225903 (1,1) 

2.25054560 (2,2) 

2.4909371 (1,2) 

2.0669861 (2,3) 



1.6510879 ± 0.619417i 1.6538681 ± 0.617998i {(. = 2) 
0.5769306 ± 0.500242i 



1.0184604 ± 0.4808314i 
0.2679495 ± 0.4999999? 



22 5.4233830 (1,1) 3.2913953 ± 0.7913153i 
5.2100924 (2,2) 1.5434373 ± 0.5190729i 
4.9289111 (3,3) 0.8664011 ± 0.500090i 



3.2918067 ± 0.7910191i (/. = 2) 
1.5531493 ± 0.5155661i (t = 4) 



5.0704584 (2,3) 1.1559456 ± 0.4978173i 
4.7922902 (3,4) 0.6426630 ± 0.4999987j 
4.5793732 (4,5) 0.2936265 it 0.5000000i 



60 14.989239 
14.957218 

14.904929 
14.834704 



(1.1) 9.3722865 ± 1.2576148i 

(2.2) 4.6566352 ± 0.6881097j 

(3.3) 3.0641980 ± 0.5397430i 

(4.4) 2.2441546 ± 0.5042061i 



9.3723056 ± 1.2575917t (C = 2) 
4.6571019 ± 0.687560i (£ = 4) 
3.0673519 ± 0.5372569i {£ = 6) 



14.931927 (2,3) 3.7935172 ± 0.1297153j 
14.871226 (3,4) 2.6121248 ± 0.4812710? 
14.793839 (4,5) 1.9631184 ± 0.4990019i 



24, 


.996095 


(1,1) 


15.745726 


± 1, 


,6103627* 


15.745730 ± 1.6103563* 


ie 


= 2) 


24, 


,984411 


(2,2) 


7.8560250 


±0, 


,84835874 


7.8561233 ± 0.8481972* 




= 4) 


24, 


,965063 


(3,3) 


5.2162758 


±0, 


,6205981j 


5.2169764 ± 0.6195799* 




= 6) 


24, 


,938297 


(4,4) 


3.8857124 


±0, 


,5339379i 


3.8883774 ± 0.5314536* 




= 8) 


24, 


,904572 


(5,5) 


3.0756442 


±0, 


,5062334i 


3.0816823 ± 0.5074819* 




= 10) 


24, 


,864499 


(6,6) 


2.5254416 


±0, 


.50068161 








24, 


,952091 


(3^4) 


4.5042779 


±0, 


.3353532* 








24, 


,922034 


(4,5) 


3.4477639 


±0, 


.4797970* 








24, 


,885222 


(5,6) 


2.7784399 


±0, 


.4976799* 









11 



Table 3: The energies and the string solutions near the breaking point for the XXX 
model. * means the exact diagonalization is not possible. 



(mi, 7712) 


N 


E exact/ J 


Ebethe/ J 


A 




91 


o.uoiyo ( ( 


o.uoiyo ( ( 


z.i 1 00 ( 00 ic u. i ( oooyoot 




22 


5 31910954 


5 31910954 


2 3727865 










2.2284619 


(2,3) 


61 


15.1840740 


15.1840740 


3.8607743 ± 0.06641234i 




62 


15.4361213 


15.4361213 


4.0193617 










3.8366348 


(3,4) 


120 


29.9665138 


29.9665138 


5.4408068 ± 0.05989996z 




121 


30.2170566 


30.2170566 


5.5399897 










5.4351240 


(4,5) 


199 


* 


49.7298474 


7.0260588 ± 0.04464858^ 




200 


* 


49.9800466 


7.1045560 










7.019490 



parity of mi + 7712. They call the string solutions as s-string for odd parity case, and 
c-string for even parity case. The string solutions are given by 

V-string(c-string): (1711,1712) = {r,r), {N — r,N — r), (0<r< — 1, (30) 



EKS-string(s-string): (mi, m2) = {s, s + 1) , {N — s — 1, N — s) , ( Sq < s < — 

(31) 

where r, s are intger and sq is minimum of mi. The distribution of sq is given by figure 
1^. For the large A^, we have 



So ~ 0.3347Ar0-4934 _ ^32) 

Therefore, we can estimate the ratio of the number of the string A^string and the total 
states A/'totai = A^(A^ - l)/2 at 

N 9 9 

i™i ~ — (33) 

A^total N nVn' ^ ' 
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Figure 1: The large behavior of the V-string sohitions for the XXX model (A = 1). 
Plot of the imaginary part of the string against \fN for {mi, 1712) = (1, 1) (plotted with 
the symbol o) and (mi, 771,2) = (2,2) (plotted with the symbol •). 




10 20 30 40 50 60 



If the string violation does not occur, the ratio is 



^^^""g ^ 1 - A (34) 
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Tabic 4: The configurations and the number of the states A^'g^j^jj^g of the complex solutions 
for the XXX model (A = 1). Here, we list for mi,m2 < N/2 cases. The other cases are 
given by taking [N — m2, N — mi). The boldface show the missing point of the string. 



N 




(mi 


m2) 




'^string/2 


8 


(1,1) 




(1,2) 




2 


1 n 










Q 
O 


20 


(1,1) 


~ (4,4) 


(1,2) 


~ (4,5) 


8 


22 


(1,1) 


~ (5,5) 


(2,3) 


~ (4,5) 


8 


60 


(1,1) 


~ (14,14) 


(2,3) 


~ (14,15) 


27 


62 


(1,1) 


~ (15,15) 


(3,4) 


~ (14,15) 


27 


120 


(1,1) 


~ (29,29) 


(3,4) 


~ (29,30) 


56 


121 


(1,1) 


~ (29,29) 


(4,5) 


~ (29,30) 


55 


199 


(1,1) 


~ (49,49) 


(4,5) 


~ (48,49) 


94 


200 


(1,1) 


(49,49) 


(5,6) 


- (48,49) 


94 


298 


(1,1) 


~ (74,74) 


(5,6) 


~ (73,74) 


143 


299 


(1,1) 


~ (73,73) 


(6,7) 


~ (73,74) 


142 


416 


(1,1) 


~ (103,103) 


(6,7) 


- (103,104) 


201 


417 


(1,1) 


~ (103,103) 


(7,8) 


~ (103,104) 


200 
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Figure 2: The large N behavior of the V-string solutions for the XXX model (A = 1). 
Plot of the real part of the string against the sites number for (mi, 7712) = (1, 1) (plotted 
with the symbol o) and (mi, 7712) = (2, 2) (plotted with the symbol •). 
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Figure 3: The large N behavior of the EKS-string solutions for the XXX model (A = 1). 
Plot of the imaginary part of the string against the sites number N for m2 = mi + 1 case. 



?31 



hC 



H 



03 

• I— I 

hC 




(6,7),(7,8), 
(5,6) 



40 60 80 100 

lattice sites number: N 



15 



Figure 4: The large N behavior of the EKS-string solutions for the XXX model (A = 1). 
Plot of the real part of the string against the sites number N for 777,2 = mi + 1 case. 



14 



12 



^ 10 
m 

H 8 



6 - 



4 - 



"1 TT" 



(1,2)..---' 



LLl. 



iiijiliiliiiiljiiililllli 



; t : t : * ! 



; > t : : i i ! : 



(3,4) 
(4,5) 

(5,6) 



20 40 60 80 100 

lattice sites number: N 



120 



140 
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in the XXX model, o means that the point is excluded. The fitting function which is 
denoted by dotted hne is given by sq ~ 7V°-^^^^. 
25 




500 1000 1500 2000 2500 3000 3500 4000 4500 5000 
lattice sites number: N 



16 



5 Conclusion 



We have presented the systematic calculations of the large N behavior in the spin 1/2 
Heisenberg XXX model. Here, we have developed a new iteration method which enables 
to evaluate all of the complex solutions. Though we have presented the calculation up to 
the site number N ~ 5000, it is possible to extend the calculation to the larger number of 

= 5000. But we believe that the essential behavior of the string solutions is clarified. 

In particular, we have studied the properties of the string solutions which are classified 
into two types, EKS-string and V-string. For the EKS-string, we find a clear evidence of 
the violation of the complex solutions, and at some number A^, they become real solutions. 
We confirm the results up to = 417. Initially, the EKS-string is a usual string solution, 
that is, the imaginary part of the EKS-string is 1/2. However, it deviates from 1/2 when 
the sites number increases, and the imaginary part of the EKS-string becomes zero. 
On the other hand, the V-string behaves in a rather normal way. We show that the real 
part and the imaginary part of the V-strings behave like A^ and V^, respectively, and 
the calculated results are consistent with the prediction by the 1/A^ expansion method. 
Finally, the complex solutions of the XXX model in the Bethe Ansatz method are classified 
by three types of solutions for given A^. The first type is given by the string solution based 
on the string hypothesis. The second is given by the EKS-string whose imaginary part is 
lower than 1/2. The third is V-string which behaves as A^" with a = 1 or 1/2. 

The EKS-string and V-string are out of the solutions based on the string hypothesis 
but not out of the Bethe Ansatz solutions. Therefore, the number of the state which is 
given by the Bethe Ansatz method is unchanged. Further, the number of the complex 
solutions is less than the prediction of the string hypothesis. But, as shown in eq.(^) 
and eq.(|3^), the difference is order Therefore, the thermodynamic quantity like as 

the entropy cannot be affected by this phenomena. 

The string hypothesis has been applied to many exactly solvable models and has 
given physically plausible results. However, even for the XXX model, some of the string 
solutions break down. Therefore, we must carefully analyze the exactly soluble model in 
the Bethe Ansatz method with the string hypothesis |T^, |T^. In the massive Thirring 



model, the bound state is calculated using the string hypothesis |18]. But the recent 



investigation shows that it is not correct |19, EG, 21 1. We will investigate for the future 
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issue how these violations of the string solutions affect the bound state problem of the 
massive Thirring or XYZ model. 
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